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We have found the exact ground state for a large class of antiferromagnetic spin- 
one chains with nearest and next-nearest neighbour interactions. The ground state 
is characterized as a matrix product of local site states and has the properties 
characteristic of the Haldane scenario. 



1 Introduction 

In recent years antiferromagnetic spin-one chains have been subject to intensive ana- 
lytical, numerical and experimental investigations. Of main interest are groundstate 
phase diagrams and critical properties with respect to variations of interaction pa- 
rameters resp. anisotropies of the systems. The study of anisotropy effects is quite 
important, as all experimental realizations indicate that quasi-one- dimensional sys- 
tems have restricted symmetries [1-7]. 

In [Q we started the investigation of a most general class of spin-one chains with 
anisotropic nearest-neighbour exchange interactions and single ion anisotropy. In a 
large parameter subspace of this model the groundstates were found in the form of 
"matrix products" . The corresponding phase diagram consists of several parts which 
are separated by transition lines of first and second order. Away from these lines 
the model has a unique ground state, an energy gap to the excited states and ex- 
ponential decay of ground state correlations. Thus, the so-called Haldane scenario 



was verified for the considered model. This scenario, conjectured in 1983 ||, [H| for 
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certain isotropic models with integral spins, is quite interesting as it points out a 
striking difference between the behaviour of isotropic integral and half-integral spin 
chains. The latter are expected to have no gap and algebraic decay of correlations. 



The first model for which the Haldane scenario was proven rigorously is the Valence- 
Bond-Solid (VBS) model — a spin-one chain with special isotropic bilinear and bi- 
quadratic nearest neighbour interactions ||TT| , [E^. In a further development the unique 
VBS ground state was cast into a matrix product of single site states and general- 
ized to anisotropic interactions [[HJ, 14j. The concept of "matrix product ground 



states" (MPGs) was then applied to most general spin-one models ||. The MPGs 
have non-trivial correlations and thus differ notably from the ground states of the 



Majumdar-Ghosh-type models [15| 



In this paper we investigate the additional effects of interactions between next-nearest 
neighbours fll6| . The anisotropies assumed in [§] for the local spin pair interactions 



will be retained. 

2 Model 

For realistic spin-one chains we assume the following symmetries: 

a) rotational invariance in the (x, y)-plane, 

b) invariance under S z — > —S z , 

c) local homogeneity of interactions, hjj + % = hj+ij. 

The most general anisotropic spin-one chain with nearest and next-nearest neigh- 
bour interactions can then be written in the following form: 

L L -y 

% = hj,j+i,j+2 = + h j+ltj+2 ) + hjj+2, (1) 

j=i i=i z 

h jij+1 = a Q A) + «! (AjBj + BjAj ) + a 2 B 2 j + a 3 A 3 - + a^Bj {1 + Bj) (2) 

+a 5 ((Sj) 2 + (^ +1 ) 2 ), 
h jJ+2 = doA] + d^AjBj + BjAj) + d 2 B) + d 3 Aj + d A Bj{\ + Bj) (3) 

+d 5 ((^) 2 +(S7 +2 ) 2 )+c, 

with real parameters ctj, dj and a constant c. We impose periodic boundary condi- 
tions. 

The nearest (next-nearest) neighbour interactions Aj and Bj {Aj and Bj) are defined 
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as 



Aj = Sj Sj~ +1 = Sj~S j+1 + Sj Sj~ +1 (transverse), (4) 

Bj = S*S* +1 (longitudinal), (5) 

Aj = Sj~Sj- +2 = Sj~Sj~ +2 + SjSj~ +2 (transverse), (6) 

Bj = S Z jS z j+2 (longitudinal), (7) 



with S ± = (S x ± iS y ) I y/2. Thus we have a model with 12 non-trivial parameters 
dj, including a scale, and an additional constant c. 



aj, 



If all the cJj's are equal to zero, we obtain the nearest neighbour model investigated 
in [|J. For this 6-parameter model it was shown that the exact groundstate could be 
found as an unique MPG in a 4-dimensional parameter sub-space. For various rea- 
sons it is interesting to know what will happen if next-nearest neighbour interactions 
are considered. The question is for which interaction parameters the MPG ansatz is 
still applicable and what is the structure of the solution manifold. 



3 Matrix product ground states 

We shall be interested in the antiferromagnetic case of the model where the ground 
state is characterized by S% otal = 0. We use the same ansatz as in || to determine 
the ground state in the form of a matrix product state. 

Denoting the Sj eigenstates with eigenvalues and ±1 by |0 >j and |± >j, we 
define at each site j a 2 x 2-matrix: 





9j = 

with non- vanishing parameters a and a. The latter will turn out to be ±1. For the 
global ground state of (ffl) we use the ansatz 

\ipo (a, cr) >= Trace (g 1 ®g 2 <g) ■ ■ ■ ®9l) , (9) 

where (g) denotes a matrix multiplication of 2 x 2-matrices with a tensor product of 
the matrix elements ||. 

By adjusting the constant c, we guarantee the condition ^j +lj+2 > (and thus 
TC > 0). We demand the MPG to be a ground state of (|T|) with eigenvalue 0: 
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TC\ipo (a, a) >= 0. This obviously is equivalent with hjj+xj+2\ipo (a, cr) >= 0. As 
is a product state, it is sufficient to demand 



(10) 



Equation (|T0|) means that the local interaction acts upon all the four entries 

of the g® 3 matrix. It can be shown that flTU| ) requires nine linear equations for the 
parameters aj and dj which guarantee the ground state property of \i/j (a, a) > under 
the condition hjj + i t j +2 > 0. As we have 12 + 2 = 14 parameters, namely the coupling 
parameters aij, dj and the two parameters a and o of the ansatz (|^), this means 
that MPGs can be constructed in a five-dimensional parameter submanifold. The 
uniqueness of the ground state (§) is achieved by a strictly four-fold degenerate lowest 
eigenvalue zero of the local interaction hjj + ij + 2 ||14j| . This condition determines the 



geometrical structure of the five dimensional solution manifold. 



4 Results 

Explicit calculations show that ( |T0"D is satisfied under the following conditions 



0) 


a = signa 3 , 




i) 


a(a + 2a 3 ) = a 3 


- «i, 


ii) 


a2 = atoCL 2 — 2 0-3 


+ 4d 3 , 


iii) 


«5 = «3 + «o(l - 


- a J — a «3 


iv) 


«0 = -«3, 




v) 


— acrolo = ct 3 — d\ 




vi) 


d 2 — (2 - a 2 )a 3 , 




vii) 


«5 = |(a 2 - 4)d 3 , 




viii) 


04 = |(a 2 — 2)d3. 





As pointed out before, the conditions (i) to (viii) define eight linear equations for the 
interaction parameters <x,- and dj whereas (0) specifies a. The most striking effect 
is that the ansatz (^) imposes strict conditions upon the parameters for the next- 
nearest neighbour interactions, equations (iv) to (viii): As soon as d-$ ^ 0, all the 
other parameters dj do not vanish either, except for special values of the parameter 
a. In particular we find: do 7^ 0. This means that as far as the next-nearest neigh- 
bour interactions are concerned biquadratic terms cannot be neglected. 



As «3 approaches zero (limit of only nearest neighbour interactions), the equations 
(iv) to (viii) become irrelevant. Besides, the d3-terms in (i), (ii) and (iii) vanish. 
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Comparing the two situations, d% = and 7^ 0, the following observation can be 
made: Starting from a model with 12 non-trivial parameters for nearest and next- 
nearest neighbour interactions, we remain with five free parameters, 03, ao, «3, «4 
and a, whereas for only nearest neighbour interactions (03 = 0) we remain with four 
free parameters, a , a 3 , a 4 and a, on the basis of a model with 6 non-trivial interac- 
tion parameters [§. 

The ground state property of (|]) can be guaranteed under the condition h > 0. 
The uniqueness of (H) is valid if all the other eigenvalues of the local interaction are 
strictly positive. The diagonalization of hjj +1 j +2 yields a set of linear, quadratic and 
cubic inequalities for the parameters aj and dj (see the appendix and [0). Those 
inequalities define a five-dimensional parameter manifold with a complex geomet- 
rical structure. In the limit dj — > this manifold turns into the four- dimensional 
sub-space defined by the conditions a > 0, a 4 > 0, a 3 7^ and o / 0. A detailed 
discussion of this solution manifold of the nearest neighbour model can be found in ||. 

For general « 3 ^ there are the following most significant effects (see also the 
appendix and [|]1|): 



We get a higher degree of freedom concerning the choice of the parameters oto an d O-a-, 
in particular they can be set equal to zero ( |A.1| ). On the other hand the ansatz (||) 
implies stricter conditions for the parameters a%, oti and thus a: The parameter |ct 3 | 
must exceed a certain finite minimum value ( |A.2| ) in contrast to the simple condition: 
\atz\ 7^ for d 3 = 0. As in || the parameter a must satisfy the condition a 7^ 0. But 
in contrast to 0, depending on the choice of the values for a and «4, a must no 
longer be chosen arbitrarily small or large (|A.2|, |A.5|). 



The ground state correlation functions can be calculated, using the transfer ma- 
trix method of 0, §]. In the thermodynamic limit L — ► 00 and for r > 2 we obtain 



the same longitudinal and transverse 2-site correlation functions with exponential 
decay as in [|J: 



1 — \a\ 



[1 - \a\y \l + \a\ 



with the longitudinal correlation length: 



e «i 



In 



1 + |a| 
1 — \a\ 



(12) 
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< SfS? >=< SfS? >= -\a\ (a + signa) 



-a 



(13) 



with the transverse correlation length 



— = In 1 1 + I a I 



(14) 



As the correlation functions depend on the overlap parameter a, the only difference 
between the cases d 3 = and a 3 ^ is given by equation (i) which describes the 
functional dependence of a on the interaction parameters. 

Note that for a —>■ there is a critical transition (diverging correlation lengths 
(0))(I1)) i 11 ^ a phase where all the spins lie in the xy-plane: \ipo >— ► |0000 . . . >. 
For a — > oo the MPG approaches the Neel states | ± =)= ± =)=... >. 



There is a finite gap to the excitations in the thermodynamic limit. This can be 
understood using the same arguments as in ||. The MPG ansatz generates a "Hal- 
dane scenario" in a natural way. The global ground state of (JT|) is composed of local 
ground states in a way which is independent of the system's size. For this reason the 
properties of the finite system, especially the energy gap, persist in the thermody- 
namic limit L —>■ oo. 



5 



Appendix 

The diagonalization of hjj+ij+2 yields the following conditions for the uniqueness of 
the MPG (g): 

a A + min(d 3 , §|a 3 |) > 0, 

«3 + ||«3| > 0, 

a + i|a 3 | > 0, (A.l) 

H > 

2a 2 d 3 + 2a A + 3|a 3 | > 0, 

a7 3 (2(a 2 -2)(2d 3 + a 4 ) + |a 3 |(a 2 + 2)) + 2a 4 |a 3 | > 0, 
2a (l + a 2 ) + 2a 2 d"3 + 3|a 3 | > 0, (A.2) 

| as | (d 3 + a ) - 2d 3 (2d 3 + a ) > 0. 

Using the following abbreviations: 



ri4 



d\ = 


2(d 3 + a 4 ), 




d 2 = 


a 4 + ||a 3 | + d 3 , 




d 3 = 


d 3 (3 — a 2 ) + |ao + «4- 


i 


d 4 = 


|a 3 | +d 3 (a 2 -2), 




d 5 = 


§(|a 3 | + a 2 a ) + d 3 (a 2 


-!)• 


d 6 = 


|"3| - as, 




d 7 = 


a + 2d 3 , 




d 8 = 


i(a + |« 3 |) + a 3 




d 9 = 






"10 = 


a 2 (a + d 3 ), 




f«o, 


n 2 = |a 3 , n 3 = 


l("3 


-«3, 


n 5 = d 3 , n 6 = 


aad 3 



(A.3) 



(A.4) 
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we further obtain: 



d 3 + d 5 

(d 3 + n 5 )(d 5 - n 5 ) - n\ 

d 3 + d 5 + d 7 + d 9 

(d 3 + d 5 )(d 7 + d 9 ) + d 3 d 5 + d 7 d 9 

+n 5 (d 3 - d 5 - n 5 ) - 2{n\ + n\ + nj) - 3nl 



> 



> 



> 



> 











(A.5) 



d 3 (d 5 + n 5 )(d 7 + dg) + d 7 d 9 (d 3 + d 5 ) 
+n|(4ni + 2n 5 ) + An 2 n 3 n 6 

-2nj(d 5 + d 9 ) - 2n\{d 3 + d 7 ) - nj(2d 3 + d 7 + 3d 9 ) 

-n 5 (d 7 + d 9 )(d 5 + n 5 ) - 2nj(d 5 + d 7 + n 5 ) + 2n^{n\ - n\) > 0. 
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